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4^ . Abstract 

I ■ We present deep J-, and i^^^-band imaging data of the MOIRCS Deep Survey (MODS), which 

2 i was carried out with Multi-Object Infrared Camera and Spectrograph (MOIRCS) mounted on the Subaru 

telescope in the GOODS-North region. The data reach 5cr total limiting magnitudes for point sources 
^ ' of J = 23.9, H = 22.8, and Ks = 22.8 (Vega magnitude) over 103 arcmin^ (wide field). In 28 arcmin^ 

of the survey area, which is ultra deep field of the MODS (deep field), the data reach the 5a depths of 
J = 24.8, H = 23.4, and Kg = 23.8. The spatial resolutions of the combined images are FWHM ~ 0.6 
arcsec and ~ 0.5 arcsec for the wide and deep fields in all bands, respectively. Combining the MODS data 
with the multi-wavelength public data taken with the HST, Spitzer, and other ground-based telescopes in 
the GOODS field, we construct a multi-wavelength photometric catalog of i^Tg-selected sources. Using the 
catalog, we present i^^-band number counts and near-infrared color distribution of the detected objects, 
CN| ■ and demonstrate some selection techniques with the NIR colors for high redshift galaxies. These data and 

. catalog are publicly available via internet. 

Key words: catalogs — galaxies :high-redshift — galaxies:photometry — infrared:galaxies — surveys 

1. Introduction break of galaxies at 1 ^ 2: ^ 4, which is important to study 

the stellar population of these galaxies and to determine 
. In recent years, several deep mult i- wavelength surveys their photometric redshift s. Compared to optical light, 
^ ' have been carried out to reveal galaxy formation and evo- observed NIR light is less sensitive to the effect of dust 
5^ , lution in the high-redshift universe. Representative exam- extinction, and the spatial resolution of NIR data is com- 
ples of such surveys are the Great Observatories Origins paratively good even in ground-based observations. These 
Deep Survey (GOODS, Giavalisco et al. 2004), the Subaru are helpful for the identification of sources detected in 
XMM-Newton Deep Survey (SXDS, K. Sekiguchi et al., other wavelengths such as X-ray, mid-IR, sub- mm and ra- 
in preparation), the Cosmic Evolution Survey (COSMOS, dio. 

Scoville et al. 2007), and the All- wavelength Extended The recent availability of wide-field NIR instruments 

Groth strip International Survey (AEGIS, Davis et al. with large- format detectors mounted on 4-8m telescopes 

2007). Deep multi- wavelength observations from radio to has allowed us to carry out wider and deeper NIR sur- 

X-ray allow us to comprehensively investigate the proper- veys. With such new instruments, we still need a rela- 

ties of stars, gas, dust, and AGN of high-redshift galaxies, tively long integration time in order to construct a rep- 

In such multi- wavelength surveys for high-redshift galax- resent at ive sample of high-redshift galaxies. For exam- 

ies, near-infrared (NIR) imaging is essential for the follow- pie, a normal L* galaxy seen in the local universe would 

ing reasons. First, the observed NIR luminosity of galaxies have i^vega ~ 23 if placed at z ~ 3 (e.g., Franx et al. 

refiects their stellar mass, which is one of the most basic 2003). Deeper data would be desirable for the investi- 

physical properties of galaxies, relatively well for galax- gation of low- mass (sub-L*) galaxies at such high red- 

ies at z ^ 3. NIR data also covers the Balmer/4000A shift. The search for star- forming galaxies at z ^ 7 also 
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Fig. 1. Survey field layout and exposure maps for J, if, and Kg bands. In the top-left panel, the GOODS-N HST/ACS region and 
the original HDF-N region are shown in magenta and green lines. Arrows represent the direction of X- and Y-axes of the mosaiced 
MOIRCS images. Grayscales in the top-right, bottom-left and bottom-right panels show the integration time of the MOIRCS J, if, 
and i^s-band imaging, respectively. 



requires extremely deep NIR data. Existing NIR surveys 
with such depth are hmited to several small field sur- 
veys. They include those with the Hubble Space Telescope 
(HST)/NICMOS such as Hubble Deep Field North (HDF- 
N; Dickinson et al. 2003; Thompson et al. 1999), HDF- 
South NICMOS field (Wilhams et al. 2000), and Hubble 
Ultra Deep Field (UDF, Thompson et al. 2005). The new 
HST instrument, WFC3 IR-channel is providing deeper 
data than those with NICMOS (e.g., Windhorst et al. 
2010), but the deep NIR observations with the HST in- 
struments are practically limited to A < 1.8/im. Although 
there are also ultra-deep NIR surveys with ground-based 
8m class telescopes such as Subaru Deep Field (Maihara et 
al. 2001) and Faint Infrared Extragalactic Survey HDF-S 
field (FIRES; Labbe et al. 2003), these surveys have small 
fields of several arcmin^. 

We have carried out a ultra-deep NIR imaging sur- 
vey, namely, MOIRCS Deep Survey (MODS) with Multi- 



Object Infrared Camera and Spectrograph (MOIRCS) 
mounted on the Subaru telescope in the GOODS-North 
region. Total integration time of ~ 124 hours were spent 
in the Ji^ i^^-band imaging observations with MOIRCS, 
and the data reach ^ 23 (Sa, Vega) over ~ 103 
arcmin^, and Kg ~ 24 over ~ 28 arcmin^. In this re- 
gion, the GOODS and the Chandra Deep Field North 
(CDF-N) surveys provided deep multi- wavelength imag- 
ing data such as optical HST/ACS images (Giavalisco et 
al. 2004), mid-IR images obtained with Spitzer/IRAC and 
MIPS (M. Dickinson et al., in preparation), and Chandra 
X-ray images (Alexander et al. 2003). Extremely deep 
radio observations with VLA, deep sub- mm/mm surveys 
with SCUBA, AzTEC and MAMBO, and extensive opti- 
cal spectroscopic surveys with the Keck telescopes have 
also been carried out (e.g., Morrison et al. 2010; Pope et 
al. 2006; Greve et al. 2008; Perera et al. 2008; Barger et 
al. 2008; Wirth et al. 2004). On the other hand, NIR 
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* the hmiting magnitude is estimated from the background fluctuation measured with a aperture diameter of 2 x FWHM 
of the PSF. The total hmiting magnitude for point sources, which is corrected for fluxes missed from the aperture, is ~ 0.3 
mag brighter for each fleld and band (see text for details). 



data which cover the GOODS-North region have reached 
only mAB ~ 22-22.5 (Capak et al. 2004; Bundy et al. 
2005), although Wang et al. (2010) recently published rel- 
atively deep i^Tg-band data which reach i^T^ ~ 22.7 (Vega). 
For J band, there had been no wide-field data which 
cover most of the GOODS-N region. Therefore NIR data 
with a comparable depth have been desirable in this field. 
In this context, with the obtained ultra-deep NIR data, 
we have investigated the number counts of Distant Red 
Galaxies at z ^ 2 (Kajisawa et al. 2006a), the cluster- 
ing properties of stellar mass-selected sample at 1 < z < 4 
(Ichikawa et al. 2007), the stellar mass dependence of 
the X-ray properties of galaxies at 2 < z < 4 (Yamada 
et al. 2009), the evolution of the galaxy stellar mass func- 
tion at 0.5 < z < 3.5 (Kajisawa et al. 2009), the correla- 
tion between the stellar mass and surface brightness for 
galaxies at 0.3 < z < 3 (Ichikawa et al. 2010), the rela- 
tion between the NIR morphology and star formation ac- 
tivity of galaxies at 0.8 < z < 1.2 (Konishi et al. 2010), 
the star formation activity as a function of stellar mass 
at 0.5 < z < 3.5 (Kajisawa et al. 2010), and the evolu- 
tion of quiescent galaxies as a function of stellar mass 
at 0.5 < z < 2.5 (Kajisawa et al. 2011). Yoshikawa et al. 
(2010) also used these NIR images with NIR spectroscopic 
data of star-forming BzK galaxies at 2: ~ 2 to investigate 
the star formation activity and stellar population of these 
galaxies. 

In this paper, we present the NIR imaging data and 
a multi-wavelength photometric catalog of iiTg-selected 
sources. Section 2 describes the observations. We give 
details of the data reduction in Section 3, and the prop- 
erties of the reduced data in Section 4. In Section 5, we 
construct the source catalog and explain the catalog en- 
tries. Using the source catalog, we present the i^g-band 
number counts and distribution of NIR colors in Section 
6. A summary is presented in Section 7. 

The Vega-referred magnitude system is used throughout 
this paper, unless stated otherwise. 



2. Observation 

The deep Ji^i^Tg-band imaging observations of the 
GOODS-North were carried out with the Multi-Object 
InfraRed Camera and Spectrograph (MOIRCS; Ichikawa 
et al. 2006; Suzuki et al. 2008) on the Subaru Telescope. 
MOIRCS consists of two 2048 x 2048 HgCdTe HAWAII-2 
detectors covering a field of view of 4 x 7 arcmin^ with 
a pixel scale of 0.117 arcsec/pixel. Our observations were 
made on 24 nights (including 8 half nights), in 2006 April- 
May, 2007 March-May, and 2008 March-May. Figure 1 
shows the field layout of the MODS in the GOODS-North 
region. Four MOIRCS pointings cover a field of ~ 105 
(7 X 15) arcmin^ (referred to as GT-1 to GT-4 field; see 
Figure 1). Nominal center positions are shown in Table 
1, and adjacent fields overlap by ~ 20 arcsec (Figure 1). 
The position angle (PA) of the field of view was fixed to 45 
degree for most of the observations. Since an engineering- 
grade chip, which is slightly less sensitive than science- 
grade one, was used on one of the channels of the MOIRCS 
(channel- 1) in 2008 March-May, however, we also used the 
PA rotated by 180 degree from the original value (i.e., PA 
= 225 degree) equally in these observations to equalize 
the depths of the fields for the both channels. For the 
iCg-band, we also used archival data of the MOIRCS ob- 
tained by other groups, namely, Wang et al. (2009) on 
2005 December and 2006 June, and Bundy et al. (2009) 
on 2006 April. 

In the observations, we used the standard circular 
dithering pattern of the MOIRCS imaging mode (1 center 
+ 8 surrounding pointings) with a diameter of ~ 12 arcsec. 
The center of the dithering pattern for each sequence was 
shifted randomly in a 10 arcsec box. Exposure times for 
each frame in the J, i^, and Kg bands are 100-260, 105- 
180, and 60-210 sec, respectively, which are split into 80- 
130, 20-40, and 20-70 sec sub-integrations, respectively, 
to avoid the saturation of the detectors (the co-adding 
mode). The total integration time for each MOIRCS 
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pointing is listed in Table 1, and the exposure maps for J, 
and Kg bands are shown in Figure 1. For the GT-1, 
3, and 4 fields, the integration times are 6.3-9.1 hours in 
J band, 2.5-4.3 hours in H band and 8.3-10.7 hours in 
Kg band, which were adjusted so that the final depth of 
the survey data is homogeneous. The GT-2 field, which 
includes the Hubble Deep Field North (HDF-N, Williams 
et al. 1996; see Figure 1), is the ultra-deep field of the 
MODS, where the total integration times are longer than 
the other fields, namely, 28.2 hours in J band, 5.7 hours 
in H band and 28.0 hours in Kg band. 

Figure 2 shows the seeing distribution of the data for 
each field. We measured FWHM of the point spread func- 
tion (PSF) for each reduced (before combined) frame us- 
ing relatively bright isolated point sources. Most of the 
data were obtained under relatively good seeing conditions 
1.0 arcsec). In excellent conditions, we preferentially 
observed the GT-2 field. 

A standard star P177-D (Leggett et al. 2006) was ob- 
served in a photometric night (1 April 2007) for flux cali- 
bration. For the i^^-band data, we used the iC-band mag- 
nitude of P177-D in Leggett et al. (2006) for the calibra- 
tion. The difference between Kg and i^-band magnitudes 
for P177-D is negligible. We obtained object frames for 
flux calibration immediately before the standard star ob- 
servation at similar airmass in all 4 pointings in JHKg 
bands. The total integration time of these object frames 
for each field and each band is 240-400 sec. 

3. Data Reduction 

The MOIRCS data were reduced basically using 
a purpose-made IRAF-based software package called 
MCSRED ^ (Tanaka et al. 2010). The reduction was per- 
formed independently for each chip, and then the com- 
bined images for the different chips (channel- 1 and 2) and 
fields (GT-1 to GT-4) were mosaiced finally. The basic 
procedure in the MCSRED is 1) making a quick object 
mask for each raw frame, 2) making a self-sky fiat frame 
by median stacking of masked raw frames, 3) flat-fielding, 
4) making a median sky frame for each flat-fielded frame 
from 6-8 frames before and after the frame, 5) perform- 
ing sky subtraction, 6) removing residual sky gradients by 
fitting each quadrant of the object-masked frame with a 
3rd-order polynomial surface and subtracting it from the 
unmasked image, 7) correcting the optical distortion, 8) 
measuring positional offsets among the frames and align- 
ing, 9) stacking the frames. Although the self-sky flat for 
NIR cameras often suffers from thermal emission from the 
camera, MOIRCS was designed so that the thermal emis- 
sion from the camera and telescope was minimum (Suzuki 
et al. 2008). We confirmed that the self-sky flat agrees 
within 5% with a dome-flat made from the "on — off" 
frame even in Kg band, where the effect of thermal emis- 
sion is expected to be maximum. We also checked the 
photometry at the different positions in the detectors us- 
ing objects in the overlap regions between the different 

^ http: / / www.naoj .org/staff/ichi/MCSRED / mcsred_e. html 



chips and fields (see below). Once the reduction sequence 
was finished, we used the combined image to make a more 
accurate object mask which includes fainer objects than 
a quick mask, and then repeated the procedures of 2)- 
9). Such "second-pass" procedure allows us to remove the 
effect of faint objects on the flat-fielding and sky subtrac- 
tion. 

In addition to the basic procedure of the MCSRED^ we 
also performed a defringing process. Each raw frame ob- 
tained in the MODS, which was carried out before the 
installation of the "fringe-free" type filter, contains a cir- 
cular fringe pattern produced by the interference of OH air 
glows. Since the amplitude and center position of the cir- 
cular pattern varied with time, the fringe pattern cannot 
be completely removed by the standard sky subtraction 
procedure. These residual patterns could systematically 
affect the source detection, photometry, or morphologi- 
cal analysis for objects in the combined images. Figure 3 
shows an example of the defringing process for a raw frame 
(channel-2, i^g-band). The circular fringe pattern can be 
seen in a frame which is flat-fielded by a dome-flat and 
is then sky-subtracted with the 3rd-order polynomial sur- 
face fitting (top- left panel of the figure). In order to reduce 
noise and highlight the fringe pattern, we masked objects 
in the image and convolved the frame with a Gaussian 
kernel with cr = 2 pixel (top- middle panel). Then we de- 
termined a center of the circular fringe pattern (e.g., yel- 
low curve in the top- middle panel), and converted the im- 
age in the Cartesian coordinate system (x, y) to that in 
the polar coordinate system (r, 6>), where the fringe pat- 
tern becomes vertical stripes (top-right panel). The con- 
verted image was smoothed along the vertical direction 
(azimuthal direction in the original Cartesian coordinate) 
with a median filter of 1 x 401 pixels to reduce noise and 
fill in the masked region (bottom-right panel). Then the 
smoothed image was converted back into the Cartesian 
coordinate system and was multiplied by the dome-flat to 
obtain a fringe frame for the original raw data (bottom- 
middle panel). Finally, the fringe frame was subtracted 
from the raw frame (bottom- left panel). 

We used this defringing process twice in the reduction 
procedure. First, we applied the defringing process to raw 
data and used the defringed frames to make a sky-flat 
frame which is not affected by the fringe pattern. With 
this sky- flat frame, we performed flat-fielding of original 
(not-defringed) raw data, and then followed the standard 
sky-subtraction procedures mentioned above. After the 
sky subtraction, we applied the same defringing process 
to the sky-subtracted frames to remove the residual of the 
time-variant fringe pattern. The defringing process after 
the sky-subtraction has the advantage of much smaller 
amplitudes of the residual patterns than in a case where 
the process is directly applied to raw data. 

Before the stacking procedure, all object frames were 
scaled to the count level of the reference frame for the 
flux calibration mentioned in the previous section, using 
relatively bright unsaturated sources. We discarded the 
frames whose count level was less than 70% of the cal- 
ibration frame, because we often failed in the sky sub- 
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Fig. 3. Example of the defringing process for a Xs-band raw frame of channel-2. top-left: a raw frame is quick flat-fielded by a 
dome-fiat and sky-subtracted, top-middle: objects in the frame is masked and then the frame is convolved with a Gaussian kernel 
to reduce noise, top-right: the frame is converted from Cartesian coordinate system (x, y) to the polar coordinate system (r, 9). 
bottom-right: the converted image is smoothed along the vertical direction with a median filter of 1 x 401 pixels, bottom-middle: 
the frame is converted back into the Cartesian coordinate system to make a fringe pattern frame, bottom-left: the fringe pattern 
frame is subtracted from the original frame. 



traction of these frames due to the bad and unstable sky 
condition. In the stacking, we used the weighted- average 
with 3(7 chpping. In order to maximize the signal-to-noise 
(S/N) ratio for the surface brightness of the combined 
image, each frame was weighted by l/cTp-^, where cTpix is 
the standard deviation in the pixel-to-pixel statistics of 
the frame whose count level was scaled to the reference 
frame. We did not weight the images with seeing size. 
Instead, we made several combined images for each field 
and each band using the data sets grouped in different see- 
ing sizes, and measured the background fluctuation of the 
combined images with apertures of various sizes. Based 
on the results, we decided to discard frames with a see- 
ing size of FWHM > 1.2 arcsec for the GT-1, 3, and 4 
flelds and > 0.8 arcsec for the GT-2 field from the stack- 
ing. For J and i^Tg-band data in the GT-2 field, we also 
provided the high-resolution combined image only from 
the data sets with a seeing size of < 0.5 arcsec to inves- 
tigate the morphology of galaxies in the NIR wavelength 
(Konishi et al. 2010). In addition to the combined image, 
we also provided the exposure map and RMS map which 
gives an estimate of the noise level at each pixel. In or- 
der to make the RMS map, we first divided the "sigma 
map" outputted from the IRAF/IMCOMBINE task by 
the square root image of the exposure map, and then nor- 
malized the resulting image so that the average of pixels 
in the RMS map corresponds to the standard deviation in 



the pixel-to-pixel statistics of the combined object image. 

In order to mosaic the combined images for the differ- 
ent chips and fields, at first, each combined image was 
rescaled to a common zero point of 26.0 AB mag (for 
all JHKs bands). We then determined the position of 
each combined image via astrometry using the GOODS 
HST/ACS version 2.0 data, whose world coordinates agree 
weh with catalogs from the SDSS, 2MASS, and deep VLA 
20cm data (Giavalisco et al. 2004; M. Giavalisco and the 
GOODS Team, 2010, in preparation). Relatively bright 
isolated stars and galaxies with early-type morphology in 
each combined image were cross- matched with those in the 
HST/ACS ^F850LP-band image. About 100 such objects 
were used for each combined image and the world coordi- 
nates were determined with a root mean square scatter of 
~ 0.03 - 0.05 arcsec, which corresponds to 0.25 - 0.45 
pixels. Based on the determined world coordinates, the 
combined images were transformed into a common grid 
of the expanded image of the GT-2 channel- 1 data. We 
performed the photometry of objects in the overlap re- 
gions between the different chips and fields to check the 
zeropoint over the detectors in each band. The magni- 
tudes at the different positions in the detectors agree well 
and no significant systematic difference is seen. The typi- 
cal measured average difference and root mean square are 
0.03 ±0.10 mag in J band, 0.02 ±0.19 mag in H band, 
and 0.01 ±0.11 mag in Kg band. In the mosaicing process. 
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the pixels in the overlap regions of the different images 
were combined with the weighted- average without sigma 
clipping. The inverse square of the RMS map for each 
combined image was used as the weight. 

4. Reduced images 

4.1. Area coverage 

The reduced J, i^^-band mosaic images have a pixel 
scale of 0.117 arcsec and a image size of 4200 x 8400 
pixels including the margin. The total exposure times 
texp for the images are shorter toward the outer edges as 
the result of the dithering technique in the observations 
mentioned above. The area in the J, and i^^-band 
images with texp > 3600 (7200) sec is 106.3 (104.3), 103.3 
(98.5), and 118.6 (114.1) arcmin^, respectively. The larger 
area coverage in Kg band reflects the contribution from 
the observations by other groups mentioned above, whose 
fleld conflgurations are different from ours. If we use the 
homogeneous depth in the central region of each fleld as 
a baseline, the area with a < 2 x cr(centrai region) the 
RMS maps is 104.6, 104.8, and 109.8 arcmin^ for the J, 
and i^g-band images. In order to make a i^^-selected 
source catalog in the next section, we use the area with 
texp > 3600 sec in the i^-band image, which ensures a < 
2 X cr(centrai region) ^11 three bauds. The area with texp > 
3600 sec in H band is 103.3 arcmin^ for the whole fleld 
and 28.2 arcmin^ for the deep GT-2 fleld. 

4-2. Image quality and PSF matching 

The PSF of the combined image in each chip and 
fleld was measured using relatively bright isolated point 
sources. The FWHMs of the PSF are listed in Table 1. 
The PSF of the J, H, and i^^-band images for the GT-1, 
3, and 4 flelds ranges 0.53 - 0.60 arcsec (0.58-0.59 arc- 
sec for the most images). The images for the GT-2 fleld, 
which were obtained preferentially in the excellent condi- 
tion, have sharper PSF with FWHMs of 0.45 - 0.49 arcsec. 

Since the different fleld images have different PSF 
FWHMs, there is PSF variance in each mosaic image. In 
order to measure colors of objects with consistent aper- 
tures across all bands, we made PSF-matched mosaic im- 
ages by convolving all images to a common PSF, namely 
that of the GT-3 i^^-band image, which has the poorest 
resolution (FWHM = 0.60 arcsec), before the flnal mosaic- 
ing. For the GT-3 J-band and GT-4 i^^-band data, we 
used the combined images provided from the frames with 
PSF FWHMs of < 1.0 arcsec for the PSF-matched mo- 
saic images in order to secure a resolution for the proper 
matching of the PSF to that of the GT-3 Kg-band image. 
Figure 4 shows a Ji^i^g-band pseudo-color image made 
from the PSF-matched mosaic images. We also made 
PSF-matched mosaic images only for the GT-2 fleld with 
a higher resolution. These images were provided from the 
combined images which were convolved to the resolution 
of the GT-2 J-band image (FWHM = 0.5 arcsec). In 
Section 5.2, we compare the growth curve of the point 
sources in the PSF-matched images among the different 
bands in order to check the accuracy of the PSF-matching. 



4-3. Background noise and limiting magnitude 

In order to estimate the background noise and derive the 
limiting magnitude, we directly measured fluxes of circular 
apertures placed at random positions on the combined im- 
ages where objects were masked and replaced with pseudo 
noise images. In order to provide the pseudo noise images, 
we averaged 10 randomly-shifted combined images (where 
objects were masked) and multiplied by the square root 
of the number of the frames used for the average in each 
pixel, assuming the Poisson statistics. We then calcu- 
lated the root mean square of these measurements, and 
repeated the measurements with various aperture diam- 
eters. Figure 5 shows the results for the different chips 
and flelds at each band. For comparison, the expected 
values from the pixel-to-pixel root mean square scaled by 
the square root of the aperture area, which is calculated 
under the assumption that the pixel-to-pixel noise is the 
Poisson distribution and uncorrelated each other, are also 
shown in the flgure (solid and dashed curves). The em- 
pirically measured background noise is larger than that 
expected from the scaled pixel-to-pixel noise, because ad- 
jacent pixels in the combined images become correlated 
through the reduction procedures such as the distortion 
correction or the registration of the frames. 

We estimated 5cr limiting magnitude for each chip and 
fleld using the background fluctuation measured with a 
aperture diameter of 2 x FWHM of the PSF. These lim- 
iting magnitudes are listed in Table 1. For GT-1, 3, and 
4 flelds, the depth of the data is relatively homogeneous 
and the limiting magnitudes are J 24.2, iJ 23.1, and 
Ks ~ 23.1. In the GT-2 fleld, the limiting magnitudes are 
J - 25.1, H - 23.7, and Ks - 24.1, which are - 1 mag- 
nitude deeper than the GT-1, 3, and 4 flelds. In order 
to derive the total limiting magnitude for point sources, 
we need to estimate the aperture correction. In Appendix 
1, we carry out simulations with the IRAF/ARTDATA 
package and flnd that the aperture magnitude with a di- 
ameter of 2 X FWHM of the PSF is fainter by - 0.28-0.29 
mag than the total magnitude for point sources. The to- 
tal limiting magnitudes for point sources corrected for the 
missed fluxes are also listed in Table 1. 

These limiting magnitudes show that the MODS data 
are one of the deepest NIR imaging data over a area of ^ 
30 to ~ 100 arcmin^. In particular, the i^^-band data in 
the GT-2 fleld is the deepest image obtained to date at the 
wavelength. For example, the ESQ GOODS-South NIR 
imaging survey with the VLT/ISAAC achieved 5<j limiting 
magnitudes of J ~ 24.0, H ~ 23.1, and Kg ~ 22.5 over 
~ 130 arcmin^, and the ultra deep i^T^-band data with a 
5<j limit of Kg ~ 23.7 was also obtained over ^6.7 arcmin^ 
in the HUDF (Retzlaff et al. 2010). The NIR imaging data 
obtained in the FIRES survey reach 5<j limits of J ~ 25.3, 
H - 24.2, and Kg - 23.8 over - 4.7 arcmin^ and J - 24.6, 
H ~ 23.7, and Kg ~ 23.1 over ~ 26.2 arcmin^ (Labbe et al. 
2003; Forster Schreiber et al. 2006). Wang et al. (2010) 
performed a wider but shallower i^T^-band imaging survey 
than the MODS in the GOODS-North region with the 
CFHT/WIRCam and achieved a bcr limit of Kg ^ 22.7 
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Fig. 4. Jif Ks-band pseudo-color image of the MODS 
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Fig. 5. Background sky fluctuation as a function of aperture diameter for each band. Different symbols represent the different 
flelds. Red and blue symbols show those for channel- 1 and 2, respectively. Dashed and solid lines show the fluctuation expected from 
the pixel-to-pixel root mean square scaled by the square root of the aperture area under the assumption of the Poisson statistics for 
the GT-1, 3, and 4 flelds and the GT-2 fleld, respectively 



over ~ 200 arcmin^. 

5. Xg-selected catalog 

We constructed multi-band photometric catalogs of Kg- 
selected sources in the MODS field, using the MODS 
data and multi-wavelength public data of the GOODS 
survey, KPNO/MOSAIC /7-band data (Capak et al. 
2004), HST/ACS B, V, i, and z-band version 2.0 data 
(Giavalisco et al., in preparation; Giavalisco et al. 2004), 
and Spitzer/IRAC 3.6, 4.5, 5.8, and 8.0 /im-band DRl 
and DR2 data (M. Dickinson et al., in preparation). Two 
catalogs are provided; one includes objects in the whole 
MODS field ("wide" catalog), and the other includes only 
those in the GT-2 field ("deep" catalog). For the deep 
catalog, the multi-band photometry was performed with 
a slightly smaller aperture (1.0 arcsec diameter) in the 
sharper PSF-matched images with FWHM of ~ 0.5 arc- 
sec than that for the wide catalog (1.2 arcsec diameter). It 
should be noted that the source detection was performed 
in the i^^-band non-convolved (not PSF-matched) image 
in both cases. We describe a method of constructing the 
catalogs and the catalog entries in the following. 

5.1. Source detection and completeness 

The source detection was performed in the i^^-band im- 
age using the SExtractor software version 2.5.0 (Bertin & 
Arnouts 1996). A detection threshold of 1.3 times the lo- 
cal background root mean square over 12 connected pix- 
els was used. In order to take account of variance of the 
background noise among different positions in the image, 
we used the RMS map as a weight image in the source 
detection procedure. The "local" -type background esti- 
mate of the SExtractor was used for the photometry. We 
adopted MAG_AUTO from the SExtractor as the total 
i^s-band magnitudes of the detected objects. In Appendix 
1, we performed simulations with the IRAF/ARTDATA 
package to estimate fluxes missed from the Kron aper- 
tures of the MAG_AUTO for point sources as a function 
of i^s-band magnitude. According to the results of the 



simulations, the missed fluxes for point sources are ex- 
pected to be ~ 2-3% dit K ^ 17-18 mag and become 
~ 10-15% near the 5<j limiting magnitudes on average. 
We also hst the magnitude corrected for these average 
missed fluxes (MAG_AUTO_COR) in the catalogs as well 
as the original MAG_AUTO. The photometric error for 
MAG_AUTO was estimated from the background fluctu- 
ation shown in Figure 5 measured with a circular aperture 
whose area is the same as the elliptical Kron aperture used 
in MAG_AUTO for each object. Since the background 
noise shown in Figure 5 is averaged over each chip, we 
also used the RMS map at the object position to take ac- 
count of the noise variation in each chip and estimate the 
background noise for the object. Since the RMS map in- 
cludes an effect of the signal from the objects, we did not 
add the photon noise from the objects to the photometric 
error separately. 

As mentioned in the previous section, we included only 
objects with tgxp > 3600 sec in the H band into the cat- 
alogs to ensure relatively uniform depth in JHKs bands. 
Total 11660 objects were detected in the selected area of 
103.3 arcmin^. In order to examine the fraction of false 
detections, we performed a simulation of the source detec- 
tion on the inverse i^^-band image with the exact same 
detection parameters. Figure 6 shows number counts of 
the objects detected in the original image and the false 
detections expected from the simulation for each chip and 
field. The 5cr limiting magnitude estimated from the back- 
ground fluctuation in the previous section is also shown 
in the figure. The expected false detection rate increases 
at faint magnitude rapidly, but it is negligible (^1 %) 
at brighter than the 5(7 limiting magnitude. Below the 
limiting magnitude (^ l-2cr of the background noise), the 
numbers of the detections in the inverse image tend to be 
larger than those in the original image. This may be be- 
cause the noise property at very faint level is affected by 
the reduction procedures such as the median sky subtrac- 
tion. 

Around very bright stars, there are also spurious sources 
caused by the image persistence in the detectors of 
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Fig. 6. Ks-band number counts of detected objects (open histogram) and expected false detections (shaded histogram) for each field 
and chip. Vertical dashed-dotted line in each panel shows the Scr limiting magnitude. Dashed line shows the detection completeness 
for the point sources. 
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MOIRCS. These spurious features are relatively faint and 
appear at ~ 8-9 arcsec from the bright stars, which cor- 
responds to the dither length between successive posi- 
tions used in the observations (Section 2). Although we 
identified 15 spurious sources around very bright stars in 
the iCg-band image and excluded them from the catalogs, 
there could remain such spurious sources in the catalogs. 
Objects near very bright stars should be handled with 
caution. 

We also estimated the detection completeness of our Kg- 
selected sample using Monte Carlo simulations with the 
IRAF/ARTDATA package. The detection completeness 
actually depends on many parameters such as magnitude, 
size, surface brightness profile, crowdedness, and so on, in 
a complex way. We here examine the completeness for the 
point source as a guide. In the simulations, a point source 
with the corresponding PSF was added at random position 
in the image. We then ran the SExtractor with the same 
detection parameters and check whether the added object 
was detected or not. For each chip and field, 200 such 
simulations were done for a given magnitude of the point 
source. The results are shown in Figure 6. The detection 
completeness for the point sources is 85-90 % level at 
the 5<j limiting magnitude. The completeness decreases 
rapidly below the limiting magnitude, and becomes ~ 50 
% at 0.5 mag below the limit. 

5.2. Color measurements 

We measured the optical-to-MIR SEDs of the de- 
tected objects, using the U (KPNO/MOSAIC), BViz 
(HST/ACS), JHKs (MOIRCS), 3.6, 4.5, 5.8, 8.0 /im 
(Spitzer/IRAC) -band images. These multi-band images 
were aligned to the i^g-band mosaic image. We con- 
volved the ACS images with the IRAF/PSFMATCH task 
to match the PSF-matched images of the MODS. Figure 
7 compares the growth curves of the point sources in the 
PSF-matched images among the different bands (ACS and 
MOIRCS bands). For the comparison, we made an aver- 
age PSF from ~ 30 relatively bright, isolated point sources 
in each PSF-matched image. The growth curves of the 
stars in the different bands agree within ± 0.02 mag at 
radii larger than 0.5 arcsec. 

For the color measurements, we basically used a fixed 
aperture diameter of 1.2 arcsec for the wide catalog and 
1.0 arcsec for the deep catalog, respectively, as mentioned 
above. For MOSAIC (/7-band) and IRAC (3.6, 4.5, 5.8, 
and 8.0 /im-band) data, whose spatial resolution is much 
poorer than that of the ACS and MOIRCS data, we first 
performed the aperture photometry for the wide (deep) 
catalog with aperture sizes of 2.28 (2.39), 2.89 (3.07), 
3.05 (3.24), 3.80 (4.02), and 4.10 (4.33) arcsec, respec- 
tively, and then applied the aperture correction by using 
the light profiles of the B- and i^^-band images convolved 
to match the resolution of /7-band and IRAC images, re- 
spectively. These aperture sizes for the U -band and IRAC 
data were selected so that the fraction of the flux sampled 
by the aperture was the same as that in the ACS and 
MOIRCS images for point sources. The SExtractor dual- 
image mode was used for the photometry. 



Before running SExtractor, we also performed the sub- 
traction of the contribution from neighboring sources for 
objects in the crowded region. For the ACS and MOIRCS 
data, we first picked up all sources in the non-convolved 
images within 12 arcsec from the i^g-selected object in 
question. We checked whether the pixel on the i^g-band 
image at the center position of each picked source belongs 
to the object in question or not in order to judge if it 
is the corresponding source. If a picked source did not 
correspond to the i^g-selected object, this source was ex- 
tracted from the non-convolved image and then was con- 
volved and subtracted from the PSF-matched (convolved) 
image. For the [/-band and IRAC data, we used the B- 
and i^5-band images for the subtraction of the contribu- 
tion from neighboring sources. In this case, all sources 
picked in the non-convolved B- or i^^-band image were 
extracted and convolved to match the [/-band or IRAC 
data at first. Then we measured a flux of the central re- 
gion of the convolved image for each picked source. Here, 
the central region is defined by the pixels belonging to the 
source in the non-convolved image. We also measured a 
flux of the corresponding (central) region in the [/-band 
or IRAC image, and scaled the flux of the convolved B- or 
i^s-band image to match that in the U -band or IRAC im- 
age. This flux scaling was done in order of the brightness 
of the sources in the [/-band or IRAC image. Once the 
scahng for the brightest source in a crowded region was 
done, the scaled image was subtracted from the image of 
the region and then the scaling for the second brightest 
source was done. After the scahng for all sources finished, 
the scaled images for sources which do not correspond to 
the object in question were subtracted from the original 
[/-band or IRAC image. 

The error for the aperture photometry was estimated 
in a similar way with that for MAG_AUTO mentioned 
above. We measured the background fluctuation for the 
convolved ACS and MOIRCS data and the MOSAIC and 
IRAC data with the aperture used in each band, and sim- 
ilarly used a RMS map for each band to estimate the 
background noise at each object position in the image. 
When the measured flux for an object is less than 2cr of 
the background noise, the upper limit is assigned. For ob- 
jects in crowded regions, we also performed the aperture 
photometry without the subtraction of the contribution 
from neighboring sources, and conservatively added a half 
of the difference between the results with and without the 
subtraction process to the photometric error. 

We included only objects which are detected above 2<j 
level in more than three bands {Kg band and other two 
bands) in the main catalogs, because it is difficult to es- 
timate the photometric redshift of those detected only in 
one or two bands. Of objects detected in the Ks band 
above 2a level, 933/10806 and 189/3983 were excluded 
from the wide and deep catalogs by this criterion. We pro- 
vide complementary catalogs of these sources separately. 
Most of these excluded objects are fainter than the 5cr 
limiting magnitude of the i^Tg-band data. 
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Fig. 7. Comparisons of the growth curves for point sources in the PSF-matched images among the different bands (BVizJHKs) 
for the wide (left) and deep (right) fields. The growth curves are normalized to mag at 2.5 arcsec radius. The growth curves in 
the top panels are plotted with vertical off'sets (by 0.1 mag) for clarity. The bottom panels show the differences from the Kg band 
in the growth curves. Vertical long-dashed lines represent the aperture radius used for the color measurements in the wide and deep 
catalogs. 



5.3. Redshifts 

Extensive spectroscopic surveys have been carried out 
in the GOODS-North region. We cross-matched our Kg- 
selected catalog with spectroscopic catalogs from the lit- 
erature. We used the catalogs by Yoshikawa et al. (2010), 
Barger et al. (2008), Reddy et al. (2006), Treu et al. 
(2005), Wirth et al. (2004), Cowie et al. (2004), Cohen 
(2001), Cohen et al. (2000), and Dawson et al. (2001). As 
mentioned above, the World Coordinate System of Kg- 
band image is based on that of the GOODS HST/ACS 
ver. 2.0 data, which agrees with these spectroscopic cat- 
alogs well. We first identified a pixel on the i^^-band im- 
age at the source position in the spectroscopic catalogs, 
and then found a i^Tg-selected object to which the pixel 
belongs. When there were different redshifts from the dif- 
ferent catalogs for a i^Tg-selected object, we chose one from 
the most recent catalog. Total 2093 i^g-selected objects 
are assigned spectroscopic redshifts. 

For all objects in the wide and deep catalogs, we also es- 
timated the photometric redshift using the UBVizJHKs^ 
3.6/im, 4.5/im, and 5.8/im-band photometry described in 
the previous section. We adopted the standard template- 
fitting technique with population synthesis models. We 
derived the photometric redshifts with three different pop- 
ulation synthesis models, namely, GALAXEV (Bruzual 
& Chariot 2003), PEGASE version 2 (Fioc & Rocca- 
Volmerange 1997), and Maraston (2005) model, and also 
used a public photometric redshift code, EAZY (Brammer 
et al. 2008) for an independent check of the photometric 
redshift. The results in all four cases are listed in the 



catalogs. 

The details of the used templates from each population 
synthesis model are described in Kajisawa et al. (2009). 
We here briefly summarize them. We adopted expo- 
nentially decaying star formation histories and assumed 
Salpeter IMF (Salpeter 1955) and Calzetti extinction law 
(Calzetti et al. 2000) for all three models. Free parame- 
ters of the models are age, the star formation timescale r, 
dust extinction, and metallicity for the GALAXEV and 
Maraston models. The model age is changed from 50 
Myr to the age of the universe at the observed redshifts, 
which is calculated for a cosmology with Ho=70 km s~^ 
Mpc-\ = 0.3 and = 0.7. For the PEGASE2, we 
used the same template sets as those used in Grazian et 
al. (2006), where eight exponentially decaying star forma- 
tion histories with various gas infalling and star formation 
timescales are assumed. The metallicity and dust extinc- 
tion were calculated in the model in a self-consistent way. 
We also added to the model templates the Lyman series 
absorption produced by the intergalactic medium follow- 
ing Madau (1995). These model SEDs were compared 
with the observed multi-band photometry and the mini- 
mum was found for each object. 

Figure 8 compares the photometric redshift with the 
spectroscopic redshift for objects with spectroscopic iden- 
tification for the cases with the different population syn- 
thesis models. The distribution of Sz (= ^hoto — ^spec) 
and the fraction of the catastrophic failure (noted as out- 
liers) with Sz/{l-\-Zspec) > 0-5 are also shown in the figure. 
The standard deviation of ^^/(l+Zspec) is similar among 
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Fig. 9. Photometric redshift distribution for Xs-selected galaxies in the wide (open histogram) and deep (shaded histogram) cata- 
logs. Four panels show the cases where different population synthesis models are used for the photometric redshift. 



the cases with the different population synthesis models 
(- 0.08-0.09 for ah redshift range and - 0.11-0.12 for ob- 
jects at ^^spec > 1-5). Although the photometric redshift 
accuracy is relatively good in any cases, the fraction of 
the catastrophic failure at z > 1.5 is slightly larger in the 
case with the Maraston model. 

Figure 9 shows the photometric redshift distribution for 
the cases with the different population synthesis models. 
In all cases, the distribution shows a primary peak at z ~ 1 
and a secondary peak at z r^^.h and a tail to high red- 
shift. The peaks in the photometric redshift distribution 
seem to correspond to the large scale structures found by 
the previous spectroscopic surveys at 2: = 0.48, 0.51, 0.56, 
0.85, and 1.02 (Cohen et al. 2000; Wirth et al. 2004). 
Although the distribution is similar among the cases with 
the different models, small systematic differences are also 
seen. In the case with EAZY, for example, the fraction of 
galaxies in the high redshift tail at z > 3 is smaller than 
the other cases, especially at z 4. On the other hand, 
the number of galaxies at z ~ 2-2.5 is slightly small in the 
case with the Maraston model. These differences seem to 
occur in the redshift range where the uncertainty of the 



photometric redshift is expected to be larger. 

54. Spitzer/MIPS 24iJ.m fluxes 

We measured 24/im fluxes for the iCg-selected objects 
using public Spitzer/MIPS 24/im data of the GOODS sur- 
vey (DR1+, M. Dickinson et al., in preparation). As de- 
scribed in Kajisawa et al. (2010), since the source con- 
fusion occurs for many objects in the MIPS 24/im data, 
we used the IRAF/DAOPHOT package (Stetson 1987) 
for the photometry in order to deal with the blended 
sources properly, following Le Floc'h et al. (2005). The 
DAOPHOT software fitted blended sources in a crowded 
region simultaneously with an PSF, which was empiri- 
cally constructed from bright isolated point sources in the 
MIPS image. The source positions in the MOIRCS Ks- 
band image are used as a prior for the centers of the fitted 
PSFs in the photometry. The photometric error was cal- 
culated from the background noise, which was estimated 
by measuring sky fluxes at random positions on the MIPS 
image. After the photometry with the DAOPHOT soft- 
ware, we also performed the aperture photometry for the 
PSF-subtracted image outputted by the DAOPHOT at 
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Fig. 10. Ks-band number counts for objects in the MODS field. Solid circles and squares show those for the wide and deep fields. 
No correction for the incompleteness or false detection is applied for the MODS data. Errorbars are based on the Poisson statistics. 
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SSDF from Minowa et al. (2005)]. 



the same source positions in order to measure the residu- 
als of the PSF fitting. For sources whose fitting residual 
was larger than the background noise, we added the fit- 
ting residual to the photometric error. In the case where 
the residual of the fitting is negligible, the 5cr limiting flux 
is ^ 20 /iJy. 

5.5. Chandra X-ray fluxes 

As described in Yamada et al. (2009), we also 
cross-matched our i^^-selected sample with the X-ray 
point source catalog of the Chandra Deep Field North 
(Alexander et al. 2003). The procedure to match the 
Chandra X-ray source positions in the catalog with the 
i^^s-selected objects was similar with that for the spectro- 
scopic catalogs mentioned above. We identified a pixel on 
the i^Tg-band image at the X-ray source position and then 
found the i^g-selected object to which the pixel belongs. 
In the case where there is no i^^-band object at the pixel, 
we extend the search pixels up to three pixel (0.35 arcsec) 
from the original pixel in order to find the counterpart 
in the i^g-band image. For most of the Chandra sources, 
the i^5-band counterparts were actually identified without 
the extended search thanks to the good positional accu- 
racy of the Chandra data. Thanks to the depth of the 
Ks-hdind data, 221 of 226 X-ray sources in the MODS 



field from the Alexander et al. (2003) 's catalog have the 
Kg-band counterparts, and only 5 sources cannot be iden- 
tified in the i^^-band image. Of these 5 sources, 4 sources 
are very faint in X-ray and are detected in only one band 
of the Chandra (full or soft or hard band) at low signifi- 
cance. Therefore these sources might be spurious sources. 
The other X-ray source is CXOHDFN J123627.5+621218, 
which is detected in all three bands (full-band flux of 3.7 
xlO~^^ ergs~^cm~'^). This source is detected marginally 
at 3.6 /im and clearly at 4.5 /im, which suggests that this 
object has a very red SED. For 221 sources with the Kg- 
band counterparts, the full, soft, and hard-band fluxes are 
listed in the Kg-selected catalogs. 

5.6. Catalog entries 

The wide and deep i^Tg-selected catalogs are publicly 
available on the MODS web site^. It should be noted 
that all magnitudes in the catalogs are in AB magnitude 
system. The conversion from the Vega system to the AB 
system in the MOIRCS JHKg bands are Jab = Jvega + 
0.915, Hab = iJvega + 1.354, and Ksab = i^^Vega + 1.834, 
respectively. The catalog entries are as follows. 

• ID — sequential identification number in order 

^ http:/ /www. astr.tohoku.ac.jp/MODS/ 
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Table 2. Xs-band number counts for the wide and deep fields. Errors of the number density are based on the Poisson statistics. 
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as outputted by SExtractor. IDs in the wide 
and deep catalogs are independently assigned (see 
ID_[deep/wide] below). 

• X_PIXEL, Y_PIXEL — pixel coordinates in the mo- 
saiced i^Tg-band image. 

• RA, Dec — right ascension and declination coor- 
dinates (J2000), which are based on the GOODS 
HST/ACS ver. 2.0 data. 

• K_AUTO — i^5-band Kron magnitude 
(MAG_AUTO from SExtractor) in AB magni- 
tude. 

• K_AUTO_COR — i^s-band Kron magnitude cor- 
rected for fluxes missed from the Kron aperture 
(assuming the surface brightness profile of point 
sources, see Section 5.1 and Appendix 1 for details). 

• SNR_AUTO — signal to noise ratio for the i^^-band 
Kron magnitude. 

• rad_AUTO — aperture radius used for the i^^-band 
Kron magnitude (2.5 x Kron radius) in pixel unit. 
The circularized value is listed. 

• rad_half — i^g-band circularized half light radius in 
pixel unit. 

• flag — extraction flags from SExtractor 

• {UBVizJHKs, 3.6, 4.5, 5.8, 8.0}_ap — aperture 
magnitude in AB magnitude in each band. The 
aperture with 1.2 (1.0) arcsec diameter is basically 
used for the wide (deep) catalog. 

• {UBVizJHKs, 3.6, 4.5, 5.8, 8.0}_err — photomet- 
ric error for the aperture magnitude in each band. 
Negative error indicates that the aperture magni- 
tude is the 2(j upper limit. 

• f_24um, err_24um — MIPS 24/im flux in /iJy and 
associated error. 

• fiux_X, flux_XSB, flux_XHB — Chandra X-ray 
fluxes in full, soft, and hard bands, respectively, in 
erg s~^ cm~^. 



• z_spec — spectroscopic redshift from the literature. 

• refer _spec — reference for the spectroscopic redshift. 
1: Yoshikawa et al. (2010), 2: Barger et al. (2008), 
3: Reddy et al. (2006), 5: Treu et al. (2005), 6: 
Wirth et al. (2004), 7: Cowie et al. (2004), 8: Cohen 
(2001), 9: Cohen et al. (2000), 10: Dawson et al. 
(2001). 

• IDref_spec — ID in the spectroscopic catalog. 

• zph_BC, zph_MR, zph_PG — photometric redshifts 
with the population synthesis models of GALAXEV, 
Maraston, and PEGASE2. 

• zph_EZ — photometric redshift outputted by 
EAZY. 

• wht_{UBViz} — weight value (inverse variance per 
pixel) at the object position in each band. 

• exp_{JHKs, 3.6, 4.5, 5.8, 8.0} — exposure time at 
the object position in seconds in each band. 

• exp_24um — exposure time at the object position 
in seconds in the MIPS 24/im image. 

• exp_{X, XSB, XHB} — effective exposure time at 
the object position in seconds in the Chandra full, 
soft, and hard-band images, respectively. 

• ID_deep, ID_wide — ID of the corresponding object 
in the other (deep or wide) catalog. 

6. Analysis 

We here present the Kg -band number counts and dis- 
tribution of NIR colors, and demonstrate some selection 
techniques with the NIR colors for high redshift galax- 
ies. In appendix 2, we also compare the i^^-band total 
magnitude and colors between the Kg band and the other 
/7-to-8.0/im bands in the MODS catalogs with those in 
other studies to check our photometry. 
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6.1. Ks-hand Number Counts 

Using the i^g-selected catalogs described above, we con- 
structed i^5-band number counts. Figure 10 and Table 2 
show the number counts for objects in the wide and deep 
fields. We limited the samples to Ks < 23 {Ks < 24) for 
objects in the wide (deep) field, where the detection com- 
pleteness is very high 90% at least for point sources) 
and the spurious detection is negligible as seen in Figure 
6. No correction for the incompleteness or the false detec- 
tion is applied. The number counts are calculated in 0.5 
mag bins and error bars are based on Poisson statistics. 
For comparison, those of other deep NIR surveys, namely, 
the UKIDSS UDS with UKIRT/WFCAM (Wilhams et 
al. 2009), the MUSYC deep NIR fields with CTIO 
4m/ISPI (Quadri et al. 2007), the GOODS-South with 
VLT/ISAAC (Retzlaff et al. 2010), the FIRES HDF-South 
with VLT/ISAAC (Labbe et al. 2003), Subaru Deep Field 
(SDF) with Subaru/CISCO (Maihara et al. 2001), and 
Subaru Super Deep Field (SSDF) with Subaru/IRCS+AO 
(Minowa et al. 2005), are also shown in the figure. Kg- 
band counts of the MODS agree well with the results of 
other surveys especially at faint magnitudes. The MODS 
number counts are slightly higher than those of other sur- 
veys at Ks ~ 18-19 and it is also seen that the number 
counts in the MODS deep field is sightly higher than those 
in the wide field. These might be due to the large scale 
structures around the HDF-N (Cohen et al. 2000; Wirth 
et al. 2004) which the GT-2 field includes as mentioned 
above. 

The logarithmic slope d{\ogN) / dm of the MODS num- 
ber counts seems to decrease with magnitude. For exam- 
ple, the best-fit slope to the number counts at 20 < i^s < 24 
is a = 0.20 ± 0.04, while the slope becomes steeper at 
Ks ^ 18. The change in the logarithmic slope of the K- 
band counts at ~ 17-18 has been reported by previous 
studies (e.g., Cristobal-Hornillos et al. 2003; Cristobal- 
Hornillos et al. 2009; Quadri et al. 2007). The slope at the 
faint magnitude of Ks > 20 is consistent with the results 
of the FIRES survey (Labbe et al. 2003; Forster Schreiber 
et al. 2006), although it is slightly flatter than those of 
the GOODS-S and SDF (Retzlaff et al. 2010; Maihara et 
al. 2001). If we divide the sample at Ks ^ 21 and esti- 
mate the slope of the number counts separately, the slope 
is a = 0.22 ± 0.06 at 18 < i^^ < 21 and a = 0.18 ± 0.05 at 
21 < i^s < 24. Although the uncertainty is rather large, 
the flattening of the slope with magnitude may continue 
even at ^ 19. 

6.2. NIR Color distribution 

Figure 11 shows J — Ks vs. Ks color-magnitude dia- 
gram for objects with Ks < 23 in the wide fleld and those 
with 23 < i^s < 24 in the deep fleld. Spectroscopically 
conflrmed stars are shown as squares in the flgure. These 
stars form a upper envelope at J — Ks ^ 0.9 in their color- 
magnitude distribution, which is consistent with the ex- 
pected locus of Galactic stars. Objects with S/N < 2 in J 
band are plotted at their 2a lower limit of the J — Ks color 
as symbols with arrows. Thanks to our deep J-band data 
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Fig. 11. J — Ks vs. Kg color- magnitude diagram for objects 
in the MODS field. Red circles show objects with Kg < 23 in 
the wide field, while blue ones show objects with 23 < Kg < 24 
in the deep field. Objects with S/N < 2 in J band are plotted 
at their 2a lower limit of the J — Kg color with arrows. Open 
squares represent spectroscopically confirmed Galactic stars. 
Horizontal dashed line shows J — Kg = 2.3, which is used as 
for the DRG selection. 

(the 2(7 limit of J ~ 25.2 for the wide field and J ~ 26.1 
for the deep field), we can do the color measurement up 
to rather red J — Ks colors. For example, we can measure 
objects with J - - 3 down to Ks - 22.2 {Ks - 23.1) 
for the wide (deep) field. The dashed line in the figure 
represents J — Ks = 2.3, which is used as the limit for the 
Distant Red Galaxies (DRGs, Franx et al. 2003). We can 
select DRGs with J-Ks> 2.3 down to Ks - 22.9 for the 
wide field and Ks ^ 23.8 for the deep field, and detected 
208 DRGs in the wide field and 81 DRGs in the deep field 
(65 are overlapping). Figure 12 shows the distribution of 
the i^5-band magnitude and redshift for the DRGs. In 
the figure, we used the spectroscopic redshift s (18 galax- 
ies) if available, otherwise the photometric redshifts with 
the GALAXEV model. The i^^-band number counts of 
the DRGs seem to turnover at Ks ^22 and do not in- 
crease with magnitude at Ks^ 22 as found in Kajisawa 
et al. (2006a). The DRGs lie mainly at 2 < z < 3.5, and 
there are tails in the redshift distribution at z ~ 0.5-2 and 
z ~ 3.5-4.5, which is consistent with the results of other 
deep surveys (e.g., Forster Schreiber et al. 2006; Quadri 
et al. 2007; Quadri et al. 2008). At low redshift {z<2), 
the number of DRGs does not increase with magnitude, 
and that of faint DRGs is relatively small. A similar trend 
is also seen at higher redshift in the magnitude range we 
investigated. On the other hand, the brightest i^^-band 
magnitude of the DRGs at each redshift gradually shifts 
to fainter values with redshift. As a result, the redshift 
distribution of the DRGs shifts to higher redshift with 
i^5-band magnitude. The DRGs at relatively low redshift 
dominate at bright magnitude {Ks ^ 19), which is also 
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Fig. 12. Redshift vs. Ks-band magnitude for DRGs in the MODS field. Red symbols show DRGs with Ks < 23 in the wide field, 
while blue symbols represent DRGs with Ks < 24 in the deep field. Upper panel shows the Xs-band number counts of the DRGs, 
and right panel shows the redshift distribution of these objects. Vertical dashed and dotted-dashed lines represent the Xg-band limit 
for the DRGs in the wide and deep fields, respectively. 
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Fig. 13. Spatial distribution of DRGs in the MODS field. 
The size of the symbols is scaled according to apparent mag- 
nitudes in the Ks-band. Red symbols show the DRGs with 
Ks < 22.9 in the wide field, while blue symbols represent those 
with 22.9 < Ks < 23.8 in the deep field. Thick solid line 
shows the survey area for which our catalog was constructed 
(texp > 3600 sec in the if -band). The original HDF-N field is 
also shown as thin line. 



consistent with the result of Consehce et al. (2007). 

Figure 13 shows the spatial distribution of the DRGs in 
the MODS field. As found by Ichikawa et al. (2007) for 
those in the GT-2 field, these DRGs seem to be strongly 
clustered. For example, several bright DRGs concentrate 
at (X~1200, Y~1900). One of these galaxies has Zgpec = 
3.661, while five of the other six DRGs show 2;phot =2.2- 
2.5 (the other one has 2;phot ^ 1-3). On the other hand, 
there are few bright DRGs in the original HDF-N field as 
shown in Kajisawa et al. (2006a). At ~ 2.5 arcmin scale 
(~ 1300 pixels), there are other similar voids of bright 
DRGs in the MODS field. The field-to-field variance of 
the number density of these galaxies seems to be large at 
the scale. 

The left panel of Figure 14 shows J — Kg vs. H — Kg 
color-color diagram for objects with Kg < 23 in the MODS 
field. We show only objects with S/N > 2 in band in 
the figure. Again, squares represent the spectroscopically 
confirmed stars. They lie mainly at 0.3 ^ J — Kg ^0.9 
and ^ H — Kg ^0.3, and there are also some red stars. 
Their distribution in the color-color diagram is consistent 
with that expected for Galactic stars, which indicates that 
the zero points of our J, and i^^-band data have been 
correctly determined. The distribution of galaxies in the 
J — Kg vs. H — Kg diagram is also consistent with those 
seen in other NIR surveys (e.g., Forster Schreiber et al. 
2006; Retzlaff et al. 2010). For an example of a color selec- 
tion technique for high-redshift galaxies in the color-color 
diagram, we show the selection criterion for Balmer break 
galaxies at z ~ 2.5 introduced by Kajisawa et al. (2006b) 
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H-Ks (Vega) Redshift 

Fig. 14. left: J - Kg vs. H -Ks two-color diagram for objects in the MODS field. Only objects with Ks < 23 and S/N > 2 in 
H band located in the wide field which are detected in H band are plotted. Objects with S/N < 2 in J band are plotted at their 
2(7 lower limit of the J — Kg color with arrows. Open squares show spectroscopically confirmed stars. Dashed-dotted line represents 
the selection criterion for Balmer break galaxies at 2; ~ 2.5 introduced by Kajisawa et al. (2006b), namely, J — Kg > 2(H — Kg) + 0.5 
and J — Kg > 1.5. right: Redshift distribution of the Balmer break galaxies selected by the JHKg color criterion. Open histogram 
shows all selected objects with Kg < 23, and shaded histogram shows those with spectroscopic identification. 



in the left panel of Figure 14 (dotted- dashed line). This 
criterion is J - Kg > 2{H - Kg) + 0.5 and J - Kg > 1.5, 
and is met while the Balmer/4000 A break of galaxies falls 
between J and H bands (i.e., 2 ^ z ^ 3), whereas most of 
the foreground and background galaxies would not meet 
such a criterion. The right panel of Figure 14 shows the 
redshift distribution of objects selected by this criterion. 
Again we used the spectroscopic redshift s if available, oth- 
erwise the photometric ones with the GALAXEV model. 
Most of these galaxies lie at 2 ^ z ^ 3.2 showing a strong 
peak at 2; ~ 2.5, although there is a relatively small con- 
tamination (~ 15 %) from low-redshift galaxies. Thus, the 
deep JHKs-hdLiid data of the MODS allow us to sample 
the SEDs around the Balmer/4000 A break of galaxies 
at 2 ^ 2: ^ 3 with high accuracy, which is important for 
investigating the stellar population of these high-redshift 
galaxies. 

7. Summary 

In this paper, we presented the deep JH Kg-hdiiid imag- 
ing data of the MODS obtained with MOIRCS on the 
Subaru telescope. The data cover an area of 103.3 arcmin^ 
in the GOODS-North region and reach the 5a total limit- 
ing magnitudes for point sources of J = 23.9, = 22.8, and 
Kg = 22.8. In 28.2 arcmin^ of the survey area, the data 
reach the 5a depths of J = 24.8, H = 23.4, Ks = 23.8. The 
World Coordinate System of the reduced images is based 
on the public HST/ACS ver. 2.0 data of the GOODS. The 
image quality of the combined images is characterized by 
a PSF with FWHM of - 0.53-0.60 arcsec for the wide 
field and ~ 0.45-0.49 arcsec for the deep field, respec- 
tively. For the color measurements, we also provided the 
PSF-matched mosaic images whose FWHM of the PSFs 



are 0.6 arcsec for the wide field and 0.5 arcsec for the deep 
field. The i^^-band detection completeness for the point 
sources is ^ 90 % at Kg ~ 23 for the wide field and at 
Kg ~ 24 for the deep field. The spurious sources are also 
negligible down to these completeness limits, although the 
false detection rate increases rapidly below the limits. 

Combining the multi- wavelength public data taken with 
the HST, Spitzer, and other ground-based telescopes in 
the GOODS field with the MODS data, we constructed 
the multi- wavelength photometric catalogs of i^g-selected 
sources. Total 9875 and 3787 objects are listed in the 
wide and deep catalogs. The catalogs also include the 
Spitzer/MIPS 24 fim fiuxes, Chandra X-ray fiuxes, and 
redshifts of the i^g-selected objects. The comparisons be- 
tween the spectroscopic and photometric redshifts suggest 
that the photometric redshift accuracy is S z / {1 -\- Zspec) ^ 
0.1 with ~ 3% outliers for galaxies with spectroscopic red- 
shifts (~ 4-12% outliers for objects with Zgpec > 1-5). 

Using the catalogs, we examined the i^^-band number 
counts and NIR color distribution. The i^T^-band num- 
ber counts in the MODS field are consistent with those 
in other general fields and show the logarithmic slope of 
d{\ogN)/dm - 0.2 at > 20. The NIR color distribution 
of spectroscopically confirmed stars indicates that the zero 
points of the MODS data have been correctly determined. 
We also demonstrated some selection techniques for high- 
redshift galaxies with the NIR colors. 

The MODS data are one of the deepest NIR imag- 
ing data over a area of ~ 30 to ~ 100 arcmin^ to date, 
especially in the Kg band. These data sample an im- 
portant part of the SEDs of galaxies at 1 ^ z ^ 4 to 
study their stellar population. The high image quality 
of the data is very useful for the deconvolution of other 
imaging data with poorer spatial resolution, while their 
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depth in the NIR is essential for the identification of ob- 
jects detected in the multi-wavelength data, which have 
been obtained by the Chandra, Spitzer, Herschel, VLA, 
and so on. These imaging data and iiTg-selected source 
catalogs are publicly available on the MODS web site 
(http: / / www.astr.tohoku.ac.jp/MODS /) . 
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Appendix 1. Simulations for estimating the aper- 
ture corrections 

In this section, we investigate how the MAG_AUTO 
from SExtractor misses fluxes as a function of magni- 
tude, using simulations with the IRAF/ARTDATA pack- 
age. Although such aperture correction could depend on 
not only magnitude but also size, surface brightness pro- 
ffie, and morphology of objects, we here examine the aper- 
ture correction for point sources as a guide. The aperture 
correction for point sources can be considered to be the 
minimum correction. The procedures are the same as in 
the measurement of the detection completeness in Section 
5.1. An artificial point source with the corresponding PSF 
was added at random position in the i^^-band image of 
each chip and field. We then ran the SExtractor with 
the same detection parameters and compared the input 
magnitude and the output MAG_AUTO. The results are 
shown in Figure 15. The MAG_AUTO is fainter by 
0.02-0.03 mag than the input total magnitude at i^T ~ 
17-18 and becomes fainter by ~ 0.15 mag near the 5<j 
limiting magnitude {K ~ 23 for the GT-1, 3, and 4 fields, 
and K ^ 24 for the GT-2 fields). The aperture correction 
seems to be slightly smaller for the GT-2 field data, which 
is ~ 1 mag deeper and has higher image quality. We fit- 
ted the difference between the input total magnitude and 
the output MAG_AUTO (i.e., the aperture correction) as 
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Fig. 15. Comparison between the input Xg-band total 
magnitude and the measured MAG_AUTO for artificial 
point sources in the simulations for each chip and field. 
Dotted-dashed line shows the best-fit second-order polyno- 
mial function. 

a function of the measured MAG_AUTO with a second- 
order polynomial function. We used the results to correct 
the MAG_AUTO and Hsted the corrected MAG_AUTO as 
MAG_AUTO_COR in our catalogs (Section 5.6). 

In the simulations, we also measured the aperture mag- 
nitude for the artificial point sources with a aperture di- 
ameter of 2 X FWHM of the PSF to estimate the aperture 
correction for the limiting magnitude measured with the 
aperture in Section 5.1. Figure 16 compares the aperture 
magnitude measured with the 2 x FWHM diameter with 
the input total magnitude. The difference between these 
two magnitudes is nearly constant over a wide range of 
the magnitude, especially at bright magnitude. We fitted 
the difference at the input Kg < 20, assuming that it is 
constant over magnitude. The resulting aperture correc- 
tion is similar among the different chips and fields and is 
~ 0.28-0.29 mag. The total limiting magnitudes for point 
sources are expected to be brighter by these values than 
those estimated from the background fiuctuation which 
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Fig. 16. Comparison between the input Kg-hand total mag- 
nitude and the aperture magnitude measured with a aperture 
diameter of 2 x FWHM of the PSF in the simulations for 
each chip and field. Dotted-dashed line represents the result 
of the linear fitting for the data at Kg < 20. 

was measured with a aperture diameter of 2 x FWHM. 
The corrected Hmiting magnitudes for point sources are 
also hsted in Table 1. 

Appendix 2. Comparison of the photometry with 
other studies in the GOODS-North 

In order to check the photometric accuracy, we here 
compare the i^Tg-band total magnitude and colors between 
the Ks band and the other [/-to-8.0/im bands in our cat- 
alogs with other studies in the GOODS-North field. The 
procedure to match sources in the catalogs of other stud- 
ies with those in the MODS catalog is the same as that 
for the spectroscopic catalogs (Section 5.3). We excluded 
objects with the deblended flag (from SExtractor) in our 
catalog from the comparisons in order to avoid the effects 
of the different spatial resolution of the images used in the 
different studies on the source detection. 

Figure 17 compares the i^^-band total magnitudes in 



the MODS catalog with those in Barger et al. (2008) 
and Bundy et al. (2009). We here used the original 
(not aperture-corrected) MAG_AUTO for direct com- 
parisons with these studies. Barger et al. (2008) ob- 
tained deep i^g-band data with CFHT/WIRCam and 
calibrated the WIRCam data using their MOIRCS data 
(Wang et al. 2009). On the other hand, Bundy et al. 
(2009) obtained relatively shallow i^^-band data with 
Subaru/MOIRCS and carried out the photometric calibra- 
tion of their MOIRCS data using the Palomar/WIRC data 
(Bundy et al. 2005). Both studies used the MAG_AUTO 
from SExtractor as the total magnitude. In the flgure, 
the i^5-band total magnitudes in the MODS catalog are 
slightly fainter on average than Barger et al. (2008) (by 
0.04 mag) and Bundy et al. (2009) (by 0.10 mag), while 
the overall agreements are relatively good in the both 
cases. Considering the zero-point uncertainty, the total 
magnitudes are consistent among the studies. 

Figure 1 8 shows comparisons oiU — Kg^ B — Kg^ V — Kg, 
i — Kg^ and z — Kg colors between Barger et al. (2008) and 
this study. It should be noted that Barger et al. (2008) 
used MAG_AUTO for the color measurements, while the 
aperture diameter of 1.2 arcsec is used for the MODS data 
(the wide catalog). Despite the difference of the aperture 
size, these colors agree well, especially in the colors be- 
tween the MOIRCS Kg band and the HST/ACS BViz 
bands. For [/-band photometry, we used the aperture cor- 
rection with the surface brightness profiles of the 5-band 
image, and the photometric error tends to be relatively 
large. The average difference of 0.12 mag in the U — Kg 
color between Barger et al. (2008) and this study is ex- 
pected from the uncertainty, and it may be due to some 
small systematic effect in the aperture correction. On 
the other hand. Figure 19 compares the B — Kg, V — Kg, 
i — Kg, and z — Kg colors between Bundy et al. 2009 and 
our catalog. We again used the aperture magnitude with 
a 1.2 arcsec diameter for the MODS data and adopted the 
aperture magnitude with a 2.0 arcsec diameter from the 
catalog of Bundy et al. (2009). An excellent agreement 
can be seen in all the colors. 

Figure 20 compares Kg -[3.6], Kg -[4.5], Kg -[5.8], 
and Kg — [8.0] between Wang et al. (2010) and this study. 
Wang et al. (2010) used a deconvolution procedure with 
a model image made from the convolved i^^-band image, 
which is similar with that used in the subtraction of the 
contribution from neighboring sources in our photometry 
(Section 5.2). On the other hand, we simply performed the 
aperture photometry after subtracting the contribution 
from neighboring sources and then applied the aperture 
correction with the surface brightness profiles of the Kg- 
band images. These colors agree relatively well within the 
uncertainty, although the colors in the MODS catalog tend 
to be bluer than those by Wang et al. (2009) at relatively 
red colors, and redder at the bluest colors, especially in 
the Kg — [3.6] and Kg — [4.5] colors. These systematic 
differences may be due to the differences in the method of 
the color measurement. 
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Fig. 17. Comparisons of the total Xs-band magnitudes (MAG_AUTO) in this study with those by Barger et al. (2008) (left) and 
by Bundy et al. (2009) (right). 
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Fig. 18. Comparisons of U — Kg, B — Kg, V — Kg, i — Kg, and z — Kg colors between this study and Barger et al. (2008). It should 
be noted that the colors from Barger et al. (2008) are measured with the Kron aperture (MAG_AUTO), while those for the MODS 
are measured with a 1.2 arcsec diameter aperture. 
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Fig. 19. Comparisons of B — , y — , i — , and 2; — Kg 
colors between this study and Bundy et al. (2009). While the 
colors in this study are measured with a 1.2 arcsec diameter 
aperture, those from Bundy et al. (2009) are measured with 
a 2.0 arcsec diameter aperture. 
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Fig. 20. Comparisons of Ks- [3.6], Xs- [4.5], Ks- [5.8], and 
Ks — [8.0] colors between Wang et al. (2010) and this study. 
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